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VERFAHREN, VORR^KJNG UND COMPUTERPROGRAMMPRODUKl^^ERZEUGUNG EINES 

DREIDIMENSIONALEN MODELLS 



BESCHREIBUN6 



Die Erfindung betrifft ein Verfahreii/ eine Vorrichtxing 
und ein Computerprograromprodukt zur Erzeugxihg eines drei- 
dimensionalen Modells fUr ein real bestehendes Objekt, 
insbesondere zur Erzeugung eines Flachen- oder Volxjmen- 
5 kOrpermodells Oder eines FE-Modells (FE: Finite Elemente) 
aus digitalisierten Daten des Objekts. 

Auf vielen technischen Gebieten besteht der dringende 
Bedarf, das Abbild eines real bestehenden Objekts in ein 

10 moglichst wirklichkeitsnahes Computermodell zu iiber- 

fUhren. Wahrend sich einfache Objekte meist noch recht • 
gut mit Hilfe von CAD-Programmen (C7UD: computerunter- 
stutztes Konstruieren) , FE-Progrcmunen oder dergleichen 
nachkonstruieren lassen, wird dies bei komplizierten 

15 Strukturen immer schwieriger. 



Besonders auffallig ist die Diskrepanz zwischen Bedarf 
und tatsachlich Machbarem in den Biowissenschaf ten, wenn 
beispielsweise die hochkomplizierte Gestalt von 

20 anatomischen Strukturen modelliert warden soil. Die 

geometrischen bzw. morphologischen Ungenauigkeiten des 
Modells koininen im besonderen MaBe dann zum Tragen, wenn 
auf der Grundlage des Modells eine FE-Analyse durch- 
gefuhrt warden soli/ um das physikalische Verhalten des 

25 Objekts zu berechnen. 

Bei der Erzeugung eines .aus Knoten und Elementen 
bestehenden FE-Modells kann prinzipiell zwischen zwei 
Methoden unterschieden werden, und zwar zwischen der 
30 direkten und der indirekten Erzeugung des FE-Modells. Bei 
der direkten Methode werden die Knotenpunkte dem FE- 
Prograinm fest vorgegeben, wahrend sich das FE-Programm 
bei der indirekten Methode die Knotenpunkte aus vorge- 
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gebenen geometrischen Elementen (z.B. aus Flachen- oder 
Volumenkorpern, Linien oder Punkten) selbst auswahlt, 

Untersucht man die bislang in den Biowissenschaf ten 
5 erstellten FE-Modelle, so fallt auf, dass die meisten der 
vorgestellten Verfahren direkte Verfahren sind, bei denen 
dem FE-Progranun die Knotenpunkte fest vorgegeben werden. 

Fig. 13 zeigt ein Beispiel far ein direkt durch 
10 Konstruktion erzeugtes FE-Modell der Gesichtsweichteile, 
wie es Motoyoshi et al • in "Finite element model of 
facial soft tissue. Effects of thickness and stiffness on 
changes following simulation of orthognatic surgery" ^ 
J Nihon Univ Sch Dent 35, Seiten 118-123 (1993) 
15 vorgestellt haben. Bei diesem Verfahren wird auf eine 
morphologisch exakte Obertragung der Objektstruktur in 
den virtuellen Raum verzichtet und es wird versucht, die 
Komplexitat der Struktur so gut wie mOglich durch 
manuelle Nachkonstruktion nachzuahmen. 

20 

FE-Modelle konnen auch direkt durch schichtweise 
Vernetzung erzeugt werden. Dabei werden die geometrischen 
Daten des Objekts durch ein schichtdiagnostisches 
Verfahren oder durch Anfertigen histologischer Schnitte 
25 gewonnen. 

Bei der punktbasierten, schichtweisen Vernetzung werden 
die gewonnen Schichtbilder digitalisiert und die Grenzen 
der interessierenden Struktur identif iziert . AnschlieBend 

30 werden in jeder Schicht auf den Grenzlinien Knotenpunkte 
definiert/ die zuerst zweidimensional in der jeweiligen 
Schichtlage und anschliefiend dreidimensional zwischen den 
einzelnen Schichtlagen vernetzt werden. Fig. 14 zeigt ein 
solches durch punktbasierte, schichtweise Vernetzung 

35 erzeugtes FE-Modell eines Zahns, wie es C. Lin et al. in 
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"Automatic finite element mesh generation for maxillary 
second premolar", Comput Methods Programs Biomed 59, 
Seiten 187-195 (1994) beschrieben haben, 

5 Die schichtweise Vernetzung kann auch voxelbasiert 

erfolgen, Dazu wird iiber jedes gewonnene Schichtbild ein 
definiertes quadratisches Gitternetz gelegt und jedem 
Quadrat, das in der Schicht einem Voxel entspricht, ein 
kubisches Element zugeordnet. Jedes Element, das in der 

10 darunter liegenden Schicht nicht zu einem gewissen 

Prozentsatz von der interessierenden Struktur bedeckt 
ist, fallt weg. Obrig bleiben Schichten aus gleich- 
formigen Elementen, die iibereinander gestapelt ein 
dreidimensionales FE-Modell aus kubischen Elementen 

15 bilden. Fig, 15 zeigt ein solches durch voxelbasierte, 
schichtweise Vernetzung entstandenes FE-Modell des 
menschlichen Schadels, wie es D. Camacho et al, in "An 
improved method for finite element mesh generation of 
geometrically complex structures with application to the 

20 skullbase", J Biomech 30, 1067-1070 (1997) beschrieben 
haben . 

Einen etwas anderen Ansatz haben P. van Zyl et al. in 
"Three-Dimensional finite element of a human mandible 

25 incorporating six osseointegrated implants for Stress 

Analysis of mandibular cantilever prostheses", Int J Oral 
Maxillofac Implants 10, Seiten 51-57 (1992) gewahlt. Das 
in Fig. 16 gezeigte FE-Modell eines Unterkiefers wurde 
direkt durch dreidimensionale Vernetzung einer aus 

30 Reflexmikroskopaufnahmen gewonnenen Punktwolke erzeugt. 
Die Punktwolke wird direkt in das FE-Programm eingegeben, 
wobei die Punkte der Punktwolke als Knotenpunkte ftir das 
FE-Modell herangezogen werden. Allerdings ist die 
automatische Vernetzung der Knotenpunkte nicht besonders 

35 zuveriassig. Die. gangigen FE-Programme haben hSufig 
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Schwierigkeiten mit der Vernetzung unterschiedlich 
dichter Punktwolken, sodass in der Kegel spezielle 
Programme zur dreidimensionalen Vernetzung der Punktwolke 
verwendet werden mUssen. Da diese Programme nicht 
5 kompatibel zu den FE-Standardprogrammen sind, entstehen 
auf der Sof twareseite Insellosungen, die nur von einem 
Fachmann bedient werden konnen. 

Indirekte Verfahren zur Erzeugung eines FE--Modells haben 
10 sich in erster Linie in den Ingenieurwissenschaf ten 

durchgesetzt • Bei den indirekten Verfahren werden dem FE- 
Programm beliebige geometrische Elemente vorgegeben, aus 
denen sich das FE-Programm selbsttatig die Lage der 
Knotenpunkte errechnet^ wobei die geometrischen Elemente 
15 lediglich die Rander und GrenzflSchen des spateren FE- 
Modells festlegen. Die geometrischen Elemente konnen 
entweder im FE-Programm selbst erstellt werden oder Uber 
eine sogenannte CAD/FEM-Kopplung als Flachen- oder 
Volumenkorper aus einem CAD-Programm import iert werden. 

20 

Um die Koordinaten der Knotenpunkte zu erstellen, gibt es 
bei der indirekten Erzeugung des FE-Modells zwei 
Mdglichkeiten: Beim sogenannten "Mapped Meshing" werden 
die Knotenpunkte so festgelegt, dass viereckige oder 
25 hexaederf5rmige Elemente gebildet werden. Im Gegensatz 
dazu werden beim "Free Meshing" dreieckige oder 
tetraederformige Elemente mit Zwischenknoten (sogenannte 
parabolische Elemente) gebildet, die sich besonders gut 
an komplexe Geometrien anpassen. 

30 

Der entscheidende Punkt beim indirekten Verfahren ist, 
dass die dem FE-Programm vorgegeben geometrischen 
Elemente in der Kegel manuell im FE-Prograram oder in 
einem CAD-Programm konstruiert werden mtlssen. In den 
35 Ingenieurwissenschaften stellt dies keineh Nachteil dar. 
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da die meisten Produkte sowieso mit Hilfe von CAD- 
Programmen entworfen werden. In den Biowissenschaf ten 
allerdings hat sich das indirekte Verfahren bis auf 
Ausnahmen^ etwa bei der Beurteilung konstruierbarer 
5 Fremdkorper wie HUf tprothesen usw., nicht durchsetzen 
kSnnen. Die Ungenauigkeiten, die entstehen, wenn ein 
reales Objekt durch direkte Konstruktion nachgeahmt wird, 
sind einfach zu groB. 

10 Den bisher einzigen Losungsweg, das Abbild eines real 
bestehenden Objekts in ein der CAD/ FEM-Kopp lung 
zugangliches Format umzuwandeln, bietet das Reverse 
Engineering. Beim Reverse Engineering wird das Objekt 
digitalisiert und werden mittels Flachenrekonstruktion 

15 CAD-Flachen erzeugt, die sich dem Grundsatz nach in ein 
FE-Programm importieren lassen. Die CAD-Flachen setzen 
sich zumeist aus frei formbaren Bezier- oder NURBS- 
Patches (NURBS: Nicht-uni forme rationale B-Splines) 
zusammen, die tiber ein Netz aus Kontrollpunkten 

20 stUckweise an die Oberf lachenf orm des Objekts angepasst 
werden. Bei den NURBS-Patches handelt es sich in der 
Regel um mindestens bikubische parametrische Fiachen- 
elemente, die die Objekt flache jeweils durch zwei poly- 
nomiale Kurven dritten Grades approximieren. Das Losen 

25 der linearen Gleichungssysteme wird mit zunehmendem Grad 
der Polynome immer aufwendiger, erlaubt aber eine 
genauere Anpassung der Flachenelemente an die Ober- 
fiachengestalt des Objekts. Aufgrund der komplizierten 
Algorithmen/ die bei der Flachenrekonstruktion zum 

30 Einsatz kommen, sind die zum Reverse Engineering 
verwendeten Sof twareprogramme bislang sehr teuer. 
AuBerdem sind die Berechnungen mit hohem Zeitaufwand 
verbunden . 
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Daneben hat sich herausgestellt, dass beim Reverse 
Engineering die Anpassung der B6zier- oder NURBS-Patches 
an die Oberf lachenform einer geometrisch komplexen 
Struktur oft sehr ungenau ist und dass dabei haufig 
5 irreguiare CAD-Flachen erzeugt werden, die sich weder mit 
CAD-Prograinmen noch mit FE-Programmen weiterverarbeiten 
lassen. 

Der Erfindung liegt die Aufgabe zugrunde/ ein Verfahren, 
10 eine Vorrichtung und ein Sof twareprogrammprodukt zur 

Verftigung zu stellen, mit denen sich mit relativ geringem 
Rechenaufwand ein genaues dreidimensionales Modell eines 
real bestehenden Objekts erzeugen lasst. 

15 Diese Aufgabe wird durch ein Verfahren gemafi Anspruch 1, 
eine Vorrichtung gemafi Anspruch 6 und ein Computer- 
prograinmprodukt gemafi Anspruch 9 gelost. 

Die Erfindung zeichnet sich dadurch aus, dass das 
20 interessierende Objekt zunachst digitalisiert wird, urn 
ein Netzmodell des Objekts zu erzeugen, das Netzmodell 
anschlieiiend in bilineare Flachenelemente zerlegt wird 
und die bilinearen Flachenelemente schlieBlich zu einem 
Flachen- oder Volumenkorpermodell vereinigt werden. 

25 

Unter dem angesprochenen Netzmodell sind Viel flachen-, 
Oberflachen- oder Polygonnetze zu verstehen, die 
typischerweise aus einer Menge endlich vieler Polygene 
bestehen, bei denen jeweils zwei Eck- oder Knotenpunkte 
30 eine Kante definieren und mehrere solcher Kanten einen 
geometrischen Korper beschreiben. Die geometrische 
Beschreibung des Korpers erfolgt in dem Netzmodell rein 
nximerisch, d.h. im Gegensatz zu einem analytischen Ansatz 
wird die geometrische Form nicht durch mathematische 
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Gleichungen definiert^ sondern allein durch die Lage und 
Dichte der Eck- bzw. Knotenpunkte. 

Da solche Netze rein numerisch beschrieben werden, lassen 
5 sie sich durch Digitalisieren des Objekts erzeugen. Das 
Digitalisieren kann dabei auf unterschiedliche Weise 
erfolgen. So kann das Objekt beispielsweise optisch oder 
beriihrend abgetastet werden, so dass sich eine die 
Objektoberf lache beschreibende Punktwolke ergibt. Aus 

10 dieser Punktwolke werden dann die Knotenpunkte fUr die 
Netzdarstellung gewonnen. Von dem zu digitalis ierenden 
Objekt k5nnen aber auch Oberflachen- oder Schnitt- 
aufnahmen angefertigt werden. Mit Hilfe dieser Aufnahmen 
werden dann die Grenzen (Boundaries) des Objekts 

15 identif iziert und werden einzelne Punkte dieser Grenzen 
wiederiam als Knotenpunkte fiir die Netzdarstellung 
verwendet . 

Bislang ging man davon aus, dass die numerischen Daten 
20 eines solchen Netzmodells nicht ohne Weiteres dazu heran- 
gezogen werden konnen, einen Flachen- oder Volumenkorper 
zu erzeugen, der durch analytische Funktionen beschrieben 
wird. So wird in der Literatur zum Beispiel davor 
gewarnt, dass in CAD-Prograinmen zwar Volumenk5rper in 
25 Oberf lachennetze konvertiert werden konnen, diese aber 
xjingekehrt nicht in Volumenkorper konvertiert werden 
konnen (vgl, Benutzerhandbuch AutoCAD 2000^, S. 719). Urn 
numerische Daten in analytische Daten xomzuwandeln, musste 
bislang der aufwendige und insbesondere bei komplexen 
30 Geometrien nicht immer erfolgreiche Weg des Reverse 
Engineering eingeschlagen werden. 

Die Erfindung schlSgt eine Brticke zwischen liumerischer 
und analytischer Beschreibung der Objektdaten, indem sie 



wo 03/102876 



PCT/EP03/05795 



- 8 - 



die niHaerischen Daten des Netzmodells in die analytischen 
Daten bilinearer Flachenelemente zerlegt- 

Unter bilinearen Flachenelementen sind FiachenstUcke zu 
5 verstehen, die jeweils durch zwei polynomiale Kurven 
ersten Grades bzw. durch zwei Strecken definiert sind. 
Die Endpunkte der Strecken ergeben sich dabei aus den 
Knotenpunkten des Netzes. Die beiden Strecken jedes 
FlSchenelements bilden zwei Kanten eines Polygonzuges, 
10 dessen ubrige Kanten sich durch Verbinden der Strecken- 

endpunkte ergeben. Jedes Flachenelement hat seine eigenen 
Kanten, die es nicht mit den angrenzenden Flachen- 
elementen teilt. 

15 Die bilinearen Flachenelemente sind vorzugsweise 

dreieckig, da sich Dreiecksf lachen besonders gut an 
komplexe Geometrien anpassen lassen. Bei einer solchen 
Dreiecksf lache fallt daher jeweils ein Endpunkt der 
beiden Strecken des Flachenelements in einem Punkt 

20 zusammen und wird die dritte Kante des Dreiecks durch die 
verbleibenden Endpunkte der Strecken gebildet. Das 
bilineare Flachenelement kann jedoch prinzipiell auch 
eine Form mit vier Kanten (z.B. ein Quadrat) einnehmen, 
bei der die Endpunkte der Strecken nicht in einem Punkt 

25 zusammenf alien. 

Sofern die bilinearen Flachenelemente durch ein CAD- 
Prograinm bearbeitet werden, sollten sie vorzugsweise in 
Form von NURBS-Patches vorliegen, da NURBS-Patches gegen- 
30 uber Rotations-, Skalierungs-, Translations- und 

Projektionsoperationen invariant sind. Dabei ist zu 
beachten, dass die Moglichkeit, NURBS-Patches hoheren 
Grades frei formen zu k5nnen, nicht genutzt wird, da die 
verwendeten Flachenelemente nur bilinear sind. 
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Um das Netzmodell in bilineare NURBS-Patches umzuwandeln, 
wird das Netz vorzugsweise in das IGES-Format (IGES: 
Initial Graphics Exchange Specification) konvertiert. Das 
IGES-Format ist ein ANSI-Standard, der ein neutrales 
5 Format ftir den Datenaustausch zwischen unterschiedlichen 
CAM- (CAM: rechnergestut zte Fertigung) und 
Computervisualisierungssystemen definiert. Die bilinearen 
Flachenelemente entsprechen dabei IGES-Elementen der 
Nummer 128, die fUr rationale B-Spline-Fiachen vorgesehen 
10 sind. 

Die einzelnen bilinearen Flachenelemente, in die das Netz 
zerlegt wurde, werden wieder vereinigt, um einen 
geschlossenen Flachenverbund oder einen geschlossen 

15 Volumenkorper zu erzeugen. Dies geschieht dadurch, dass 
die gegenuberliegenden Kanten von zwei benachbarten 
Flachenelementen zusammengehef tet werden. Es werden also 
die zuvor getrennten Kanten der Flachenelemente zu einer 
gemeinsamen Kante zusammengef asst, sodass das eine 

20 Flachenelement unmittelbar in das andere Flachenelement 
tlbergeht. Da samtliche Flachenelemente eben sind, gehen 
sie nicht stetig ineinander liber und es entsteht ein 
facettierter Flachenverbund. 

25 Dieser facettierte Flachenverbund stellt ein Flachen- 

modell des digitalisierten Objekts dar. Wenn der Flachen- 
verbund der Oberflache ein endliches Volumen quasi 
wasserdicht umschlielit, so entsteht ein Volumenkorper- 
modell des digitalisierten Korpers. 

30 

Das Fiachen- oder Volxamenkorpermodell lasst sich 
problemlos durch CAD/ FEM-Kopp lung in ein FE-Programm 
importieren und zu einem FE-Modell vernetzen, mit dessen 
Hilfe physikalische Berechnungen durchgefUhrt werden 
35 kOnnen. 
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Die Erfindung hat den Vorteil, dass sie ftir den Obergang 
zwischen der numerischen und analytischen Beschreibung 
der Objektdaten Kurvengleichungen ersten Grades heran- 
5 zieht* Im Vergleich zuiri Reverse Engineering^ bei dem 

Kurvengleichungen dritten oder hoheren Grades verwendet 
werden, reduziert sich dadurch die Komplexitat des zu 
losenden Gleichungssystems erheblich, Auch dann, wenn von 
einem sehr feinen Netz ausgegangen wird/ um ein besonders 

10 genaues dreidimensionales Modell des interessierenden 

Objekts zu erzeugen, ist die durch das weniger komplexe 
Gleichungssystem erzielte Zeitersparnis so grofi, dass 
sich der Rechenaufwand trotz der hohen Genauigkeit des 
Modells insgesamt verringert. Mit Hilfe der Erfindung 

15 lasst sich also mit relativ geringem Rechenaufwand ein 
genaues dreidimensionales Fiachen-, Volumenkorper- oder 
FE -Modell eines real bestehenden Objekts erzeugen. 

Die Erfindung kann sowohl in Form eines Verfahrens als 
20 auch in Form einer Vorrichtung oder eines Software- 
programmprodukts umgesetzt werden. 

Bei der Vorrichtung ist zu beachten, dass neben einer 
Datenverarbeitungseinrichtung, die die Daten- 

25 verarbeitungsschritte Einlesen des Netzmodells, Zerlegen 
des Netzmodells in bilineare Fiachenelemente, Vereinigen 
der bilinearen Flachenelemente zu einem Flachen- oder 
Volumenkorpermodell und gegebenenf alls Erstellen eines 
FE-Modells aus dem Flachen- oder Volumenkorpermodell 

30 ausftihrt, auch eine Digitalisierungseinrichtung 

vorgesehen ist, mit der sich das Netzmodell des Objekts 
erzeugen lasst. Unter eine solche Digitalisierungs- 
einrichtung fallen samtliche bildgebenden Gerate wie 
Fotoapparate und RSntgengerate, die zweidimensionale 

35 analoge oder digitale Bilder erzeugen, aus denen in 
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Kombination mit einer Bildverarbeitung ein dreidimensio- 
nales Netz gewonnen werden kann. Unter eine solche 
Digitalisierungseinrichtung fallen aber auch optische und 
beruhrende Abtastgerate, die die Oberflache des Objekts 
5 abtasten, um elne dreidimensionale Punktwolke zu 

erzeugen, aus der dann in Kombination mit einer Bild- 
verarbeitung ein Netz gewonnen wird- In diesem Zusammen- 
hang ist es ohne Belang, ob die Digitalisierungs- 
einrichtung und die Datenverarbeitungseinrichtung 
10 r^umlich getrennt sind oder ob die Digitalisierungs- 
einrichtung die Datenverarbeitungseinrichtung zur 
Ausftihrung der Bildverarbeitungsschritte nutzt. 

Es liegt auf der Hand, dass die Erfindung anstatt durch 
15 die Vorrichtung auch durch ein Computerprogrammprodukt 
realisiert werden kann, das die angesprochenen Daten- 
verarbeitungsschritte anhand von Sof twareroutinen 
abarbeitet, wenn es auf einem Computer lauft. Das 
Computerprogrammprodukt kann auf einem Datentrager 
20 gespeichert sein oder direkt in den Arbeitsspeicher des 
Computers geladen werden. 

Weitere Aufgaben, Vorteile und Merkmale der Erfindung 
ergeben sich aus der folgenden ausfuhrlichen Beschreibung 
25 von AusfUhrungs- und Vergleichsbeispielen. Dabei wird auf 
die beigefvigten Zeichnungen Bezug genommen, die Folgendes 
zeigen: 

Fig. 1 ein durch optische Abtastung digitalisiertes 
30 Abbild eines Probandengesichts in Form einer Punktwolke; 

Fig. 2 die Punktwolke von Fig. 1, nachdem sie ausgedunnt 
wurde; 



35 



Fig. 3 ein aus der Punktwolke von Fig. 2 erzeugtes Netz; 
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Fig. 4 ein Ausschnitt aus dem in Fig, 3 gezeigten Netz; 

Fig* 5a bis Fig. 5d eine erlauternde Darstellung, wie aus 
5 drei Polygonen des in Fig. 4 gezeigten Ausschnitts drei 
bilineare Flachenelemente gebildet und wieder vereinigt 
we r den; 

Fig. 6 ein aus dem in Fig. 3 gezeigten Netz erzeugtes 
10 Flachenmodell; 

Fig. 7 ein aus dem in Fig. 6 gezeigten Flachenmodell 
erzeugtes FE-Modell der Gesichtsweichteile; 

15 Fig. 8 eine digitalisierte Rontgenschichtaufnahme eines 
menschlichen SchSdels in Hohe des Unterkief ers, in denen 
die Schadelgrenzen markiert und mit Punkten versehen 
sind; 

20 Fig. 9 ein aus mehreren Schichtaufnahmen erzeugtes 

digitales Abbild eines Schadels in Form einer Punktwolke, 
nachdem sie homogenisiert wurde; 

Fig. 10 ein aus der Punktwolke von Fig. 9 erzeugtes Netz; 

25 

Fig. 11 ein aus dem in Fig. 10 gezeigten Netz erzeugtes 
Volumenkarpermodell des Schadels; 

Fig. 12a ein Netzmodell des menschlichen Ohres und 
30 Fig. 12b ein aus diesem Netzmodell durch Reverse 
Engineering erzeugtes Modell aus NURBS -FlSchen; 

Fig. 13 ein durch Nachkonstruktion erzeugtes FE-Modell 
von Gesichtsweichteilen gemaB Stand der Technik; 
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Fig. 14 ein durch punktbasierte, schichtweise Vernetzung 
erzeugtes FE-Modell eines Zahns gemaB Stand der Technik; 

Fig. 15 ein durch voxelbasierte/ schichtweise Vernetziing 
5 erzeugtes FE-Modell eines menschlichen Schadels gemaB 
Stand der Technik; und 

Fig. 16 ein durch Vernetzung einer Punktwolke erzeugtes 
FE-Modell eines Unterkiefers gemaB Stand der Technik. 

10 

Anhand der Figuren 1-7 wird nun ein erstes Ausfuhrungs- 
beispiel beschrieben, wie sich aus einem digitalisierten 
Abbild eines Probandengesichts ein dreidimensionales 
Flachenmodell und ein FE-Modell der menschlichen 
15 Gesichtsweichteile erzeugen lassen. 

Das Gesicht eines Probanden wurde mit Hilfe eines licht- 
codierten Triangulationsverf ahrens (TRICOLITE™ der Firma 
Steinbichler) digitalisiert . Dabei wurde durch einen LCD- 

20 Projektor eine Serie von Streif enmustern auf das Gesicht 
geworfen, die durch zwei CCD-Kameras aus unter- 
schiedlichen Blickwinkeln erfasst wurde. Der komplette 
Messvorgang dauerte etwa zwei Sekunden. Durch 
geometrische Auswertung (Triangulationsprinzip) wurde 

25 daraus ein dreidimensionales Abbild der Gesichts- 

oberfiache in Form einer Punktwolke erzielt. Nahere 
Einzelheiten zu diesem lichtcodierten Triangulations- 
verf ahren werden der Dissertation (citation in progress) 
von Holberg "iErfassung von Gesichtsoberf lachen durch 

30 ein lichtcodiertes Triangulationsverf ahren", Ludwig- 

Maximilians-Universitat Mtinchen (2002) zu entnehmen sein. 

Die erhaltene Punktwolke wurde anschlieBend gefiltert, urn 
eine bestimmte Aufldsung zu erzielen und redundante Daten 
35 einzusparen. Und zwar wurden die Bildpunkte/ deren Lage 
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gegeniiber den benachbarten Bildpunkten nur wenig abwich^ 
geloscht^ wahrend die Bildpunkte, deren Lage gegentlber 
den benachbarten Bildpunkten starker abwich, beibehalten 
wurden, Dadurch ergab sich die in Fig. 2 gezeigte 
5 ausgedunnte Punktwolke, in der die Bildpunkte umso 
dichter liegen, je starker sich die Topologie der 
Gesichtsoberf lache andert. 

Die ausgedunnte Punktwolke wurde anschlieliend in das 
10 Bildbearbeitungsprogramm Rapid Form™ (INUS Technology, 

Inc.) importiert und zu einem Polygonnetz aus dreieckigen 
Polygonen vernetzt. Um ein durchgehendes^ sauberes 
Polygonnetz zu erhalten, wurde das Polygonnetz von 
Lochern und von kreuzenden, redundanten und nichtmannig- 
15 faltigen Flachen befreit. Im vorliegenden Fall wurde die 
Moglichkeit, das Polygonnetz zu homogenisieren oder die 
Anzahl der Polygonf lachen zu erhohen oder zu reduzieren, 
nicht genutzt. Es entstand das in Fig. 3 gezeigte 
lochfreie und gereinigte Polygonnetz, das im DXF- Format 
20 zwischengespeichert wurde. 

Das im DXF~Format abgespeicherte Polygonnetz wurde 
anschlieliend in das Programm PolyTrans'^ (Okino Computer 
Graphics) importiert, um das Polygonnetz, ohne welter 

25 verandert zu werden, im neutralen IGES-Format 

abzuspeichern, das den Datenaustausch zwischen unter- 
schiedlichen CAD-, CAM- und Computervisualisierungs- 
systemen erlaubt. Durch den Export als IGES-Datei wurde 
das Polygonnetz in bilineare NURBS-Patches mit der 

30 Elementnummer 128 zerlegt. Die entstandene IGES-Datei 
wurde anschlieBend in das CAD-Programm Mechanical 
Desktop'^ (Autodesk Inc.) eingelesen, wobei das Abbild des 
Probandengesichts nunmehr in Form von einzelnen Flachen- 
elementen (bilinearen NURBS-Patches) vorlag, die jeweils 

35 einem Polygon des ursprtlnglichen Polygonnetzes 
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entsprachen. Mit der Funktion "Flachen zusanunenhef ten" 
(in anderen Programmen auch als "ZusammenfUgen" Oder 
"Stitchen" bezeichnet) wurden die einzelnen Flachen- 
elemente dann wieder zu einem Flachenverbund vereinigt/ 
sodass das in Fig. 6 gezeigte Flachenmodell des 
Probandengesichts entstand. Dieses Flachenmodell wurde 
dann uber die CAD/FEM-Schnittstelle AMACISOUT des 
Programms im SAT-Format export iert, um die Geometrie des 
Flachenmodells der CAD/FEM-Kopplung zuganglich zu machen. 

Die einzelnen Verarbeitungsschritte, denen die Bilddaten 
beim Zerlegen des Polygonnetzes in die einzelnen Flachen- 
elemente und beim Vereinigen der Flachenelemente 
unterlagen, werden nun ausfUhrlicher anhand von Fig. 4 
und den Figuren 5a bis 5d erlautert. 

Fig. 4 zeigt einen Ausschnitt des in Fig. 3 gezeigten 
Polygonnetzes, der dem rechten Wangenbereich des 
Patienten entnommen wurde. In diesem Bereich werden nur 
die drei fett markierten Polygone betrachtet, deren Eck- 
bzw, Knotenpunkte hervorgehoben sind. Diese drei Polygone 
sind in Fig. 5a ohne ihre Umgebung dargestellt. 

Die in Fig. 5a dargestellten Polygone wurden durch die 
Konvertierung ins IGES- Format in die drei in Fig. 5b 
gezeigten bilinearen NURBS-Patches umgewandelt. Unter 
NURBS-Patches sind wie gesagt Flachenelemente zu 
verstehen, die jeweils durch zwei nicht-uni forme 
rationale B-Splines definiert sind, d.h. durch zwei frei 
formbare polynomiale Kurven. Da die Polygone, von denen 
ausgegangen wurde, eben waren, sind die beiden B-Splines, 
die jedes Flachenelement definieren, ebenfalls nicht 
gekriimmt und entsprechen somit Kurven erstes Grades. Bei 
der Erfindung wird also jedes NURBS-Patch nicht mehr 
niomerisch durch die Eckpunkte des jeweiligen Flachen- 
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elements, sondern analytisch durch zwei Kurven ersten 
Grades bzw. durch zwei Strecken beschrieben. Da die 
ursprtinglichen Polygene dreieckig waren^ sind auch die 
entstandenen NURBS-Flachen dreieckig. Das bedeutet, dass 
5 jeweils ein Endpunkt der beiden Strecken des NURBS- 
Patches in einem Punkt zusammenfallt und die beiden 
Strecken zwei Kanten des dreieckigen NURBS-Patches 
bilden/ wahrend sich die dritte Kante durch Verbinden der 
beiden verbliebenen Streckenendpunkte ergibt. Da jeder 
10 NURBS-Patch durch ein eignes Paar linearer B-Splines 
definiert ist, hat jedes NURBS-Patch seine eigenen 
Kanten, die es nicht mit den angrenzenden Flachen- 
elementen teilt. 



15 Die einzelnen NURBS-Patches wurden dann durch die 

Zusammenhef tfunktion wie in Fig. 5c und Fig. 5d gezeigt 
jeweils mit den benachbarten NURBS-Patches vereinigt. Das 
Zusammenhef ten erfUllt dabei zwei Aufgaben: Zum einen 
werden mit dieser Funktion zwei oder mehr zusammen- 

20 hangende Fiachen zusammengehef tet, urn einen Flachen- 

verbund zu erstellen, und zum anderen werden Fehler in 
der Geometrie oder Topologie, die wahrend der 
Konvertierung aufgrund von unterschiedlichen internen 
Toleranzen und Berechnungsf ehlern auftreten kOnnen, 

25 korrigiert. Durch das Zusammenheften entsteht also ein 
durchgangiger FlStchenverbund, der als Fiachenmodell 
herangezogen werden kann. 

Das in Fig. 6 gezeigte, tiber die CAD/FEM-Schnittstelle 
30 export ierte Fiachenmodell des Probandengesichts wurde in 
das FE-Programm Design Space™ von T^sys, Inc. importiert. 
Das Flachenmodell wurde dabei wie ein gekrtimmtes Flachen- 
tragwerk behandelt, das sich nach der Schalentheorie 
verhait. Nach Zuweisung einer einheitlichen Dicke wurde 
35 das Konstrukt zu dem in Fig. 7 gezeigten dreidimensiona- 
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len FE-Modell der Gesichtsweichteile vernetzt. Die 
Vernetzung erfolgte ohne Schwierigkeiten, da fehlerhafte 
Flachen bereits wahrend der Reinigung des Polygonnetzes 
beseitigt warden, sodass es zu keinen Oberschneidungen 
5 kam. Das FE-Modell wies im GroBen und Ganzen die gleiche 
hohe Auflosung wie das Polygonnetz und das im CAD- 
Prograinm erstellte Flachenmodell auf . Durch Festlegung 
entsprechender Lagerungen und Lasten lieBen sich mit 
diesem FE-Modell hochauf losend Verformungen, Spannungen 
10 und Dehnungen in den Gesichtsweichteilen berechnen* Diese 
Berechnungen kOnnen zum Bei spiel ftir die Planung 
kosmetischer Operationen genutzt werden. 



Anhand der Figuren 8-11 wird nun ein zweites Ausftihrungs- 
15 beispiel beschrieben, wie sich aus digitalisierten 

Schichtbildern eines menschlichen Schadels ein drei- 
dimensionales Volximenkorpermodell erzeugen iSsst. 

Von dem Schadel eines Probanden wurde mit Hilfe eines 
20 Computertomografieverfahrens ein Satz von zweiundvierzig 
digitalen Rontgenschichtbildern gewonnen* Mit Hilfe des 
Bildbearbeitungsprogramms 3D-Doctor™ (Able Software) 
wurden in den Schichtaufnahmen die Grenzen des Schadel- 
knochens identif iziert und die markierten Grenzlinien mit 
25 mehreren Punkten versehen. Fig. 8 zeigt ein Beispiel 
eines solchen mit Grenzlinien und Punkten versehenen 
Schichtbilds im Bereich des Unterkieferknochens • 



Die in den einzelnen Schichtbildern markierten Punkte 
30 wurden dann zu einer dreidimensionalen Punktwolke 

zusammengefasst, die redundanten Daten durch Filtern mit 
einer bestimmten AuflOsung entfernt und die verbleibende 
Punktwolke homogenisiert . Daraus ergab sich die in Fig. 9 
gezeigte ausgedtinnte Punktwolke, in der die Bildpunkte 
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umso dichter liegen, je starker sich die Topologie der 
Schadeloberf lache andert. 

Die ausgediinnte Punktwolke wurde anschliefiend zu einem 
5 Polygonnetz aus dreieckigen Polygonen vernetzt und das 
Polygonnetz von kreuzenden, redundanten und nichtmannig- 
faltigen Flachen gereinigt. Eventuell entstandene Locher 
wurden wieder geschlossen. Das sich ergebende Polygonnetz 
ist in Fig. 10 gezeigt. 

10 

Das Polygonnetz wurde dann ahnlich wie im ersten 
Ausftihrungsbeispiel in das Programia PolyTrans™ (Okino 
Computer Graphics) importiert und in bilineare NURBS- 
Patches mit der Elementnummer 128 zerlegt, die 

15 anschlieiJend wieder unter Korrektur der auf getretenen 

Geometrie- oder Topologief ehler zu einem Flachenverbund 
zusammengehef tet wurden. Da bereits das Polygonnetz frei 
von Lochern war und der Flachenverbund daher der 
durchgehenden Oberflache eines in sich geschlossenen 

20 festen Korpers (Solid) entsprach, wurde beim Zusammen- 
heften automatisch ein Voliamenkorpermodell erzeugt* Das 
fertige Volumenkorpermodell ist in Fig. 11 gezeigt und 
hat im GroBen und Ganzen die gleiche hohe Aufl5sung wie 
das ursprUngliche Polygonnetz, sodass es fUr ein hoch- 

25 auflOsendes FE-Modell geeignet ist. Das entstandene FE- 
Modell konnte zur Simulation der Auswirkungen von 
Gewalteinwirkungen auf den Schadel herangezogen werden. 

Die obigen AusfUhrungsbeispiele zeigen, dass sich mit der 
30 Erfindung hochgenaue dreidimensionale Modelle real 
bestehender Objekte erzeugen lassen. 



Ein Vergleich der Erfindung mit dem herkOmmlichen 
Verfahren des Reverse Engineering soil aufierdem den 
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Nachweis bringen, dass sich mit der Erfindung auch der 
Rechenaufwand senken lasst. 

Hierzu wurde von einem aus 338 6 dreieckigen Polygonen 
5 bestehenden Netzmodell des menschlichen Ohres ausgegan- 
gen. Das in Fig. 12a gezeigte Netz wurde zum einen nach 
den MaBgaben der Erfindung und zum anderen mit Hilfe der 
Reverse-Engineering-Funktion des Programms Rapid Form"* 
(INUS Technology, Inc.) in NURBS-Patches der IGES- 

10 Elementenummer 128 umgewandelt. In beiden Fallen warden 
die 338 6 Polygene in ebenso viele NURBS-Patches 
umgewandelt* Im Fall des Reverse Engineering wurde 
auJierdem ein Modell mit einer geringeren Auflosung von 
232 NURBS-Patches angefertigt, das in Fig. 12b gezeigt 

15 ist. Die Berechnungszeit und der Speicherbedarf fUr die 
einzelnen Modelle lassen sich Tabelle 1 entnehmen. 



Tabelle 1 





NURBS 
Anzahl 


Rechen- 
zeit 


Speicherbedarf 
( xinkortprimiert ) 


Speicherbedarf 
( konpr imiert ) 


Erfindung 


3386 


2,9 s 


1,2 m 


0,4 MB 


Reverse 
Engineering 


3386 


1983 s 


29,6 MB 


16,1 MB 


Reverse 
Engineering 


232 


98 s 


2,0 m 


1,1 MB 



20 

Tabelle 1 ist zu entnehmen, dass das erf indungsgemafie 
Verfahren verglichen mit dem Reverse Engineering bei 
gleicher Auflosung (d.h. bei gleicher Anzahl an NURBS- 
Patches) deutlich weniger Rechenzeit und Speicherbedarf 
25 erfordert. Wenn beim Reverse Engineering die Anzahl an 
NURBS-Patches gesenkt wird, lassen sich zwar die Rechen- 
zeit und der Speicherbedarf reduzieren, doch werden 
dadurch geometrische Ungenauigkeiten in Kauf genommen* 



wo 03/102876 



«20- 



PCT/EP03/05795 



Wie der Vergleich der Figuren 12a und 12b ergibt, erzeugt 
das Reverse Engineering zwar sehr glatte OberflSchen- 
tlbergange, doch kann es insbesondere im Randbereich des 
Modells zu fehlerhaften Darstellungen kommen (siehe Pfeil 
5 in Fig. 12b) . Die Ungenauigkeiten im Randbereich sind 
lamso auffalliger, je komplexer die nachziabildende 
Struktur ist und je geringer beim Reverse Engineering die 
Aufiesung (Anzahl der erzeugten NURBS-Patches) ist. 
Vergleichbare Darstellungsf ehler treten bei der Erfindung 
10 nicht auf. 



Dartiber hinaus stellte sich heraus, dass die durch das 
Reverse Engineering erzeugten NURBS-Patches insbesondere 
dann, wenn das Patch-Layout automatisch erstellt wurde, 

15 zu irregularen, fehlerhaften und vertwisteten CAD-Flachen 
fUhrten, die nicht (z.B. tiber CAD/ FEM-Kopp lung) weiter- 
verarbeitet werden konnten. Da die manuelle Erstellung 
des Patchy-Layouts zeitaufwendig, f ehlertrachtig und in 
der Regel nur vom Faclimann zu bewerkstelligen ist, bietet 

20 das Reverse Engineering keine praxisgerechte Alternative, 
um aus digitalisierten anatomischen Strukturen wie dem in 
Fig. 12a gezeigten Ohr FE-Modelle zu erzeugen. 

Wie sich aus der obigen Beschreibung ergibt, eignet sich 
25 die Erfindung nicht nur als hochgenaues, Zeit xind 

Speicherplatz einsparendes Alternatiwerf ahren fUr das 
Reverse Engineering, sondern erschlieBt die Erfindung 
auch neue Anwendungsbereiche bei der Model lierung 
hochkomplexer Objekte, etwa in den Biowisserischaf ten, die 
30 dem Reverse Engineering und anderen herk5mmlichen 
Verf ahren bislang nicht zuganglich waren. 



wo 03/102876 




PCT/EP03/05795 



PATENTANSPRUCHE 

1. Verfahren zur Erzeugung eines dreidiitiensionalen 
Modells fur ein real bestehendes Objekt, mit den 
Schritten: 

Digitalisieren des Objekts, uiti ein Netzmodell des 
5 Objekts zu erzeugen; 

Zerlegen des Netzmodells in bilineare Flachen- 
elemente; und 

Vereinigen der bilinearen Fiachenelemente zu einem 
Fiachen- Oder Volumenk5rpermodell des Objekts. 

10 

2. Verfahren nach Anspruch 1, bei dem die bilinearen 
Flachenelemente dreieckig sind, 

3. Verfahren nach Anspruch 1 oder 2, bei dem die 
15 bilinearen Flachenelemente NURBS-Patches sind. 

4. Verfahren nach Anspruch 1, mit dem Schritt Erstellen 
eines Finite-Elemente-Modells aus dem Flachen- oder 
Volumenkorpermodell . 



wo 03/102876 PCT/EP03/05795 

22 



5. Verfahren nach Anspruch 1, bei dem das Netzmodell 
aus Punktwolkendaten des Objekts gewonnen wird, 

5 6. Vorrichtung zur Erzeugung eines dreidimensionalen 
Modells ftir ein real bestehendes Objekt, mit: 

einer Digitalisierungseinrichtung, mit der sich ein 
Netzmodell des Objekts erzeugen lasst; und 

einer Datenverarbeitungseinrichtung zur Durchfuhrung 
10 folgender Datenverarbeitungsschritte : Einlesen des 
Netzmodells; Zerlegen des Netzmodells in bilineare 
Fl^chenelemente; und Vereinigen der bilinearen Flachen- 
elemente zu einem Flachen- oder Volumenkorpermodell des 
Objekts. 

15 

7. Vorrichtung nach Anspruch 6, bei der die Daten- 
verarbeitungseinrichtung aus dem Flachen- oder Volumen- 
korpermodell ein Finite-Elemente-Modell des Objekts 
erstellt. 

20 

8. Vorrichtung nach Anspruch 6 oder 1, bei der die 
Datenverarbeitungsschritte in der Datenverarbeitungs- 
einrichtung anhand von Sof twareroutinen abgearbeitet 
werden. 

25 

9. Computerprogrammprodukt zur Erzeugung eines 
dreidimensionalen Modells fUr ein real bestehendes 
Objekt/ das anhand von Sof twareroutinen folgende Daten- 
verarbeitungsschritte abarbeitet, wenn es auf einem 

30 Computer lauft: Einlesen eines Netzmodells des Objekts; 
Zerlegen des Netzmodells in bilineare Fiachenelemente; 
und Vereinigen der bilinearen Fiachenelemente zu einem 
Flachen- oder Volumenkorpermodell des Objekts. 
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10. Computerprogrammprodukt nach Anspruch 9, das die 
bilinearen Flachenelemente erzeugt, indem es das 
Netzmodell des Objekts durch Konvertierung ins IGES- 
Fonuat in bilineare NURBS-Patches zerlegt. 

11. Computerprogrammprodukt nach Anspruch 10/ bei dem 
das ins IGES-Format konvertierte Netzmodell 
ausschliefllich IGES-Fiachenelemente der Nummer 128 
enthalt . 

12. Computerprogrammprodukt nach Anspruch 9, das aus dem 
Flachen- oder Volumenkorpermodell durch CAD/FEM-Kopplung 
ein Finite-^Elemente-Modell des Objekts erstellt. 
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Fig. 12a 




Fig. 12b 
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